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Abstract 

The effect that the current quark mass Mo may result in nonzero-ness of chiral 
condensates is systematically reexamined and analyzed in a two-flavor Nambu-Jona- 
Lasinio (NJL) model which simulates Quantum Chromodynamics (QCD) at temper¬ 
ature T and finite quark chemical potential fi without and with electrical neutrality 
(EN) condition and at any T and without EN condition. By means of a quantitative 
investigation of the order parameter m indicating the chiral symmetry breaking in the 
model’s ground states, it is shown that a nonzero Mq is bound to lead to nonzero quark- 
antiquark condensates throughout chiral phase transitions within the frame of the NJL 
model, no matter whether the order parameter m varies discontinuously or continu¬ 
ously. In fact, a complete disappearance of the quark-antiquark condensates are proven 
to demand the non-physical and unrealistic conditions /i > or 3> Mq if T = 0 

and finite, or T —>■ oo if /r < -|- Mq , where A is the 3D momentum cut of the loop 

integrals, the largest physical mass scale in the NJL model. Theoretically these results 
show that when Mq is included, besides the explicit chiral symmetry breaking indicated 
by Mo , different from the chiral limit case, we never have a complete restoration of 
dynamical (spontaneous) chiral symmetry breaking, including after a first order chiral 
phase transition at low T and high fj.. In physical reality, it is argued that these results 
play a decisive role in the known phase diagram of the model. It is the nonzero-ness of 
the quark-antiquark condensates that leads to the appearance of a critical end point in 
the first order phase transition line and the crossover behavior at high T and/or high fi 
cases, rather than a possible tricritical point and a second order phase transition line. 
They also provide the basic reason for that one must consider the interplay between the 
chiral and diquark condensates in the research on color superconductor at zero T and 
high case. The whole discussions make us learn how a source term of the Lagrangian 
(at present i.e. the current quark mass term) can greatly affect dynamical behavior of 
a physical system. 


1 INTRODUCTION 

By means of the Nambu-Jona-Lasinio (NJL) model [I] which simulates Quantum Chro¬ 
modynamics (QCD) , one has made extensive research on chiral phase transitions at fi¬ 
nite temperature T and finite quark chemical potential fj, (corresponding to quark matter 
density) [H |3]. It is found that in the chiral limit where quarks have zero current masses, 
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the dynamical (spontaneous) chiral symmetry breaking in vacuum induced by the quark- 
antiquark condensates, or say, chiral condensates (qq) [U lU El |6] will be restored at high 
T and/or high /i, and this restoration will always be accompanied with the condensates 
(qq) = 0, no matter whether the symmetry restoring phase transition is first or second order 

[ZlIH]- 

On the other hand , when there is a nonzero current quark mass Mq, in the vacuum we 
will have both the explicit chiral symmetry breaking indicated by Mq and the dynamical 
chiral symmetry breaking induced by the condensates {qq) ^ 0, thus the order parameter 
m indicating chiral symmetry breaking will include the contributions from both Mq and 
the condensates {qq). Usually one views m as a whole order parameter indicating chiral 
symmetry breaking and does not distinguish it into Mq and the contribution from the con¬ 
densates {qq). However, in the process of chiral phase transitions at high T and/or high 
/r, seeing that Mq is a fixed parameter, what can actually be changed is merely the {qq) 
sector. Therefore, the chiral phase transitions practically only depend on variation of the 
condensates {qq), but such variation must presuppose existence of the current quark mass 
Mq. In the Lagrangian of the model, the current quark mass Mq corresponds a source term 
whose existence, in general, will inevitably affect on dynamical behavior of the model, in 
present case, i.e. the change of the condensates {qq). In view of this, it is certainly an 
significant topic of the change of the chiral condensates {qq) and its effect on chiral phase 
transitions when a current quark mass Mq exists. To focus our attention on this topic and 
to compare the obtained results with the ones in the chiral limit case where Mq = 0 will be 
able to make us acquire a deeper theoretical insight of chiral phase transitions of the NJL 
model. 

As stated above, in the chiral limit, the condensates {qq) could be equal to zeros at some 
high T and/or high /i, hence the dynamical chiral symmetry breaking will be restored com¬ 
pletely. Now one can ask if a similar result could appear when the current quark mass Mq 
exists, i.e. at some high T and/or high /x, the condensates {qq) disappear totally thus the 
dynamical (spontaneous) sector of the chiral symmetry breaking will be restored completely. 
If the answer is no, then what physical effects this will actually bring about? 

The above questions motivate the research of the present chapter. Although there have 
been many work about chiral phase transitions in a NJL model with a current quark mass, 
it seems still to lack a systematical and deep examination of the problem from the point of 
view considering the relation between the current quark mass and the chiral condensates, 
or say, the one between the explicit and spontaneous chiral symmetry breaking. The dis¬ 
cussions in present chapter will just focus to this respect. It should be indicated that in the 
discussions we will revisit the process of chiral phase transitions at different conditions and 
this will naturally involve, at least qualitatively and partly, some known results. However, it 
is essential that we will analyzes all the the results from the above specific point of view and 
anticipate to get a deeper understanding about the mechanism of chiral phase transitions 
of the NJL model. 

In the following discussions , we will use the term ’’complete restoration of the dynam¬ 
ical chiral symmetry breaking” to express the condensates {qq) = 0. However, it should 
be remembered that in the case with the current quark mass Mq, different from the chiral 
limit case , the complete restoration of dynamical chiral symmetry breaking does not mean 
chiral symmetry restoration, because after all we always have the explicit chiral symmetry 
breaking indicated by Mq. 

We will explore the above topic by means of a two-flavor NJL model which simulates 
QCD in the mean field approximation. The Lagrangian of the model can be expressed by 

£ = Mo)q +Gsiiqq)'^ + {qn5rq)'^] (1) 
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with the quark Dirac fields q in the SUf{2) doublet and the SUc{3) triplets, i.e. 


q = 


Ui 

di 


i = r,g,b (three colors), 


T = (Ti,r 2 ,r 3 ) are the Pauli matrices, Gs is the four-fermion coupling constants and Mq 
represents the common current mass of the u and d quarks. The chiral SUfL{2) ® SUfR{2) 
flavor symmetry of the Lagrangian (1) will be broken not only explicitly by the current quark 
massMg, but also assumedly in the vacuum spontaneously by the scalar quark-antiquark 
condensates (qq) formed through the four-fermion interactions Gs{qq)'^- The constituent 
quark mass i.e. the order parameter indicating chiral symmetry breaking will be m = 
Mq — 2Gs{qq)- It is seen from this definition of m that the chiral (qq) = 0 will mean that 
m = Mq. The emphasis of research will be put on the effect of the current quark mass 
Mq on changes of the chiral condensates {qq) and the relevant physical results, including 
the interesting problem that if the dynamical chiral symmetry breaking could be restored 
completely at a high T and/or a high g when Mq exists. 

In Sectl2]and[31 we will describe the model’s chiral phase transitions in T = 0 and high 
fj, case respectively without and with electric neutrality (EN) condition and compare the 
obtained results. By means of a quantitative analysis of the locations of the least value 
points of the effective potential, we will focus on the changes of the chiral condensates in 
the phase transitions and the comparison of the results with the ones in chiral limit case. 
We will also find out the conditions in which the chiral condensates could be reduced to 
zeros and indicate the non-physical feature of the conditions. In Sect01 a general analysis 
of chiral phase transitions of the model will be conducted in any T and /i cases without EN 
condition. We will indicate the decisive role of nonzero-ness of the chiral condensates in the 
phase diagram and give a general demonstration of the nonzero-ness of the chiral condensates 
within the frame of the NJL model. Finally, in Sectl5]we come to our conclusion. 


2 ZERO T AND HIGH /i PHASE TRANSITIONS 
WITHOUT ELECTRICAL NEUTRALITY 


In this section, we will revisit the chiral phase transitions of the model at T = 0 and finite /r 
without the EN condition and track the variations of the order parameter indicating chiral 
symmetry breaking as g increases. 

When r = 0, in the mean field approximation, the effective potential of the model can 
be expressed by [5] 


V{m,n) 

(m - Mo)^ 


- 12 


/ \/l^ + M§ + 0it.i-E){g-E) , E = ^jP + m2 


4Gs 

(m — Mq)^ 3 


4Gs 


+9{g - m) 


2'k'^ 


A(A^ -I- - — I K\J A^ -I- w? + In 


A -1- VA^ -I- 


— (m —>■ Mo) 


^ \/— 4m^) -|- [ ^\/-b In 

o Z \ 


/i +fE — w? 


( 2 ) 


where m is the constituent quark mass i.e. the order parameter indicating chiral symmetry 
breaking, ^ is the common chemical potential of the two flavor quarks and A is 
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the 3D momentum cutoff of the loop integrals. We will focus on the least value points of 
the effective potential V(m, which correspond to the ground states of the model. In all 
the following calculations, A and the four-fermion coupling constant Gs will be fixed by the 
formula of the tt decay constant 




d^p 6{K^ — ff) 

(27r)‘* (p2 _ 77j2 _|_ 

A 


3m^ A -|- VA^ -I- m? 

^ I m 


VA^ -I- 


Gell-Mann Oakes Renner relation [101 


2 Mo{qq) 
m^ = -— 


(3) 

(4) 


and the vacuum form of the following gap equation from the extreme value condition 
dV{m,p)/dm = 0 expressed by 


m = Mq H- 


A\/h? -I- m? — rn? In 


A -f y/ A? + w? 


m 


-9{fi -m)(- m2 - In ^ ^ - 

\ m 


(5) 


By inputting the experimental value m,n.=139.57 MeV m, the phenomenological value 
/^=92.4 MeV [H] and taking the current quark mass 


Mo = 5.55 MeV, 


( 6 ) 


we will obtain from Eqs.(3),(4) and the equation (5) with p = 0 (vacuum) 

A = 636.944 MeV, Gs = 5.2866 x 10"® MeV"^. (7) 

In the meantime we also get m = = 322.39 MeV which is the only minimum point of 

V(m, p) at p — 0, since it is easy to verify that the second derivation of V(m, p) over m 


d'^V 

dm? 


1 3 r A^ o 2 A A -1- VA2 -f m2 A \ 

2Gs TT^ y/Js? + rn^ + w ^ n ^ v'A^ -|- j 

+«(p - m) - 3 In 


( 8 ) 


satisfies 


d^V 


dm'^ 


= 35779.3 MeV^ > 0. 


m—m^ 


This means (both explicit and spontaneous) chiral symmetry breaking in vacuum. When 
we increase p, the extreme value points of V(m,p) will be determined by Eqs.(5) and (8). 
The results have been shown in Table 1. Based on them we may discuss the chiral phase 
transitions of the model at T = 0 and a finite p. 
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Table 1 Variations of the extreme value points of V(m, fj.) as rising of the quark chemical 
potential fi. Denotations {V2)i and V respectively represent the values of d^Vfdm^ and V 
at the extreme value points (* = 1,2,3). 


(MeV) 

0 

300 

337 

339.22 

340 

350 

500 

636.97 

to)' 

(MeV) 

(V2)i 

(MeV^) 

322.39 

322.39 

308.56 

301.47 

297.85 




> 0 

> 0 

> 0 

> 0 

> 0 




Vi/10® 

(MeV^) 



-6.901 

-6.915 

-6.921 




TO2 

(MeV) 



174.11 

211.71 

223.10 




(H 2)2 

(MeV^) 



< 0 

< 0 

< 0 




V 2 /IO® 

(MeV^) 



-6.763 

-6.858 

- 6.886 




n 

m 3 

(MeV) 



124.17 

102.37 

97.62 

65.78 

11.12 

5.55 

(Mg) 

(^ 2)3 

(MeV^) 



> 0 

> 0 

> 0 

> 0 

> 0 

> 0 

V 3 /IO® 

(MeV'^) 



-6.768 

-6.915 

-6.969 

-7.736 

-31.66 

-83.38 


First we note that, if we increase /x but keep /x < m)* = 322.39 MeV, for instance, 
/X = 300 MeV, then V(m,/i) will have the only minimum point still at m = m)*, as shown 
in the left-third column. This indicates that the chiral symmetry breaking in vacuum will 
be maintained. However, if ^ goes up to above m®, then the situation will be changed. 
In the left-fourth column of Table 1 we give an example of /x = 337 MeV > rrii. In this 
case, the gap equation (5) will have three solutions corresponding to two minimum points 
m = = 308.56 MeV, = 124.17 MeV and a maximal point m = = 174.11 MeV. 

Obviously , since V(to)') < V{m^), to)' = 308.56 MeV will still be the least value point 
of V (to, /x), though it has been less than to® = 322.39 MeV in vacuum. It is found that as /x 
increases further, V{mi ) — V{m^) will go up, and finally it will occur that H(to)') > V{m'^) 
so that the least value point of V (to, transfers from to)' to TOg discontinuously, as shown 
in the /x = 340 MeV column of Table 1 and this indicates that a first order phase transition 
has happened. The critical chemical potential /Xd may be determined by the condition 

V (to)'‘=' , /Xd) = H (to^'^ , /Xd) 

where to)'°' and TOg”"^ are respectively the two solutions of Eq.(5) with d^Vldm? > 0. From 
this we obtain /Xd = 339.22 MeV with = 301.47 MeV and TOg”' = 102.37 MeV, 
as shown in the left-fifth column of Table 1. It is emphasized that after the the Hrst 
order phase transition, we have the least value point TOg Mg and this implies that the 
quark-antiquark condensates {qq) remain to have quite large contribution to to, hence the 
dynamical sector of chiral symmetry breaking, to quite large extent, has still not been 
restored. As a comparison, we note that in the chiral limit (Mg = 0), we always have the 
V’ corresponding second minimal point TOg = 0 [3, hence the present TOg ^ Mg can only 
be attributed to the current quark mass effect. Of course, due to existence of Mg, we could 
at most have rxig = Mg after a hrst order phase transition. Assume that is the case, then 
the situation will be a little similar to the chiral limit case, i.e. through a hrst order phase 
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transition, the sector of the dynamical chiral symmetry breaking will be restored completely, 
though the explicit chiral symmetry breaking indicated by Mq remains. Then the least value 
point m 3 = Mo will no longer change as a further increase of /r. The phase transition will 
cease at the first order critical value /r = /Xd. However, the result Mq has obviously 

negated the above assumption. 

In fact, as /i continues to increases, the least value point of V (m, fj.) will smoothly 
decrease and at above some value of /i, it will become the only minimum point of V (m, fj.) 
left, as shown in the right-third column. A direct question is whether it may be expected 
that as p, grows up further to some critical value iic 2 , we will have t Mq, thus the 

dynamical chiral symmetry breaking will be restored completely. To answer this question, 
we may take m = Mq in Eq.(5) and obtain the equation to determine fic2 

+ (9) 

A* + V M - ^0 

whose solution is 

fj, = = + = 636.97 MeV, (10) 

in view of Eqs.( 6 ) and (7), as is shown in the right-first column. However, it is noted that the 
resulting ’’second order” critical chemical potential ^c 2 has exceeded the 3D momentum cut 
off A which should be considered as the reasonable largest mass scale in this NJL model and 
is nonphysical. The above result shows that when the current quark mass Mq exists, within 
the frame of the NJL model where ^ < A, the limit m —>■ Mg or (qq) = 0 is not realizable. 
This fact implies, on the one hand, that the dynamical chiral symmetry breaking can not 
be restored completely; and on the other hand, that after a first order phase transition one 
will get only a smooth crossover behavior of the order parameter m which, however, can 
never arrive at the limit value Mg. 


3 ZERO T AND HIGH /i PHASE TRANSITIONS 
WITH ELECTRICAL NEUTRALITY 

In the above discussions we did not consider the electrical neutrality (EN) condition of the 
quark matter. In this section we will extend our discussions to the case with EN condition 
and examine if similar conclusion can be reached. It is noted that the quark matter with 
EN is a more realistic case where electron must be included in the chemical equilibrium of 
weak decays of the quarks. In the electric neutrality case, the effective potential V (m, /x, /Xg) 
in the mean field approximation can be expressed by [5] 
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V[m,n,iJ.e) 

_ (m - Mof 
4Gs 
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2{E -JlP+ M^) + e{pu - E)ipu - E) 


+^(Md “ E){p4 — E) 


127r2’ 


E = ^p^ + w? 


{m — MqY 3 

4Gs 47r2 


<1 2A(A^ + — rr? ^AV^A^ + vn? + rr? In 


A + VA^ + rn?- 


-{m Mo) 


+ (^OiPu - m) {fJ,l - 4:m^) + ^ (/iuV/iS - rn 


+771,2 In ■ 


where p = —dVIdn is the quark chemical potential corresponding to the total quark number 
density n, pe is the chemical potential of electron and 

2 1 

Pu — p ^ pe: Pd — E ~ Pe — pu 4“ +ej (i^) 

are respectively the chemical potential of the u, d quarks. The last equality of Eq.(12) is 
usually referred as beta equilibrium [9] 

For deriving the electrical neutrality condition, it is noted that the electrical charge 
density in the two-flavor quark matter with electrons is 

2 1 


with nu, Ud and Ue denoting respectively the number density of the m, d quark and electron. 
From it we may obtain 

_ _ dV duQ 

due dnq due 

Hence the FN condition will become ng = —dV/dpg = dV/dpe = 0 and has the following 
explicit expression 


dV 

dpe 


5^ 


/ 2/ie.2 2 

(+-—) -m 


1,12 




Me N 2 


1 3/2 


= 0 . 


(13) 


Owing to Fq.(13), for a given /i, we will have pe, = peim). Hence, in view of the EN 
condition dV/dpe = 0, the extreme value condition of the effective potential V(m,p,pe) 
now becomes 

dV dpe dV ^ 

dm dm dpe dm dm ’ 
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then from the first equality of Eq.(ll), it is easy to obtain the gap equation dV/dm = 0 
with the explicit form 


m = Mq H-5- m ■ 


2 ( A \/In 


A + V 


- m) - rri^ In ^ ™ ^ Md) | ■ (14) 


Owing to dV/dm = dVidm in the EN condition, the second derivative of V{m, fj,, He) over 
m may be expressed by 

fV_ _ d ^ d /dV\ dne 

dm? dm \ dm J d^e V ^m J dm 

d'^V d'^V dne 

dm^ d^edm dm ' 

From the EN condition dV/d^e = 0 we may have 

-f 

dm \ 


dV\ _ d^V d^V dne 
d^e) dmdfie 9/ig dm 


It leads to 


dne 


d^V ^d^V 


dm dmd^e ^ d^l 

Hence the second derivative of V(jn, fi, fie) over m becomes 

fV_ _ d^V / d'^V V d'^V 

dm^ dm^ \dmdfie) dfi^ 
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H-^ < 3m^ In 


A + VA^ + 


A 


A3 


m 


•\/A^ + j VA^ + 


2Gs 7r2 

^ 9{fiu -m) (fiu y/fil-m"^ - 3to^ In ^ 'l + (m« 9d) 


+ 


29{iiy, - m)y/fil - m2 - 9{fid - m)^Jfi^ - m? 


2 ^ 


W{fiu - m)fiu^/fil - to2 + 9{fid - m)fid^/fi^-m"^ + 3fil 


. 


(15) 


Now for a given /r, the extreme value points of V{m, fi, fie) will be determined by the si¬ 
multaneous equations (13) and (14), and except for this, the whole discussions of the chiral 
phase transitions under EN condition will be parallel to the ones made in SectlJJ In Table 
2 we list the variations of the extreme value points of V{m, fi, fie) for some selected values 
of ^ in a successively increasing order. 










































Table 2 Variations of the extreme value points of V{m, fj,, fj,e) as rising of /i under 
electrical neutrality condition. Denotations /iei, {dV2)i and Vi (i = 1,2,3) represent 
respectively the values of /Xe, cPV/dm? and V at the extreme value points mf. 


(MeV) 

0 

300 

345 

345.6 

(Mci) 

346 

350 

500 

656.87 

(Mc2) 

n 

rrii 

(MeV) 

322.39 

322.39 

283.69 

277.81 

272.27 




Mel 

(MeV) 

0 

0 

23.89 

26.18 

28.23 




idV2)i 

(MeV^) 

> 0 

> 0 

> 0 

> 0 

> 0 




Vi/IO® 

(MeV*) 



-6.958 

-6.967 

-6.974 




u 

m 2 

(MeV) 



189.12 

209.32 

221.93 




Me2 

(MeV) 



52.31 

47.32 

43.97 




{dV2)2 

(MeV^) 



< 0 

< 0 

< 0 




V 2 /IO® 

(MeV*) 



-6.928 

-6.953 

-6.968 




u 

m 3 

(MeV) 



143.58 

130.74 

125.3 

96.37 

12.3 

5.55 

(Mo) 

Me3 

(MeV) 

{dV2)3 

(MeV^) 



62.28 

64.53 

65.54 

70.65 

109.41 

143.81 



> 0 

> 0 

> 0 

> 0 

> 0 

> 0 

V 3 /IO® 

(MeV*) 



-6.931 

-6.967 

-6.992 





The left-second and left-third column of Table 2 indicate that, and as has been actually 
checked, the chiral symmetry breaking in vacuum will be maintained up to /i < m)* = 
322.39 MeV. If ^ = 345 MeV > to)* in vacuum, then the effective potential 
will show two minimal points m)* and mg and a maximal point m^. It may be seen that 
in this case Vi < V3 and this merely implies that the least value point rrii of V{m, fj,, ^e) 
in vacuum smoothly changes to a lower value and no phase transition happens. However, 
if /X = 346 MeV, we will have Vi > V 3 . This indicates that the ground state has been 
transferred discontinuously to the minimal point (mg = 125.3 MeV, fies = 65.54 MeV) 
and a first phase transition must have happened. The critical chemical potential /x = 
should be determined by the equation 


V(m^, /X, /Xei) = V(m^, /X, ^es), 

where (m(^,/Xei) and are respectively the two minimal points of V(m, /i,/Xe). The 

result /X = = 345.6 MeV is listed in the left-fifth column of Table 2. 

At this point, the smaller minimal value point mg ^ = 130.74MeV ^ Mq, this again 
implies that after the first order phase transition the chiral condensates {qq) have not been 
reduced to zeros, thus the dynamical chiral symmetry breaking has not restored completely. 
It is seen from Table 2 that as /x continue to go up, for example, for /x > 350 MeV, the effec¬ 
tive potential V (m, /x, /Xe) will have the only minimal point mg left whose values will decrease 
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gradually. Finally mathematically = Mq can be arrived at ^ = 656.87 MeV and 

Me = Me 3 = 143.81 MeV which are solutions of Eqs.(13) and (14) with m replaced by Mq, 
as shown in the left-last column of Table 2. The smooth change of the quark mass to 
its current mass Mq indicates that complete restoration of the dynamical chiral symmetry 
breaking, if it is realizable, seems to be a ’’second order phase transition”. However, it is 
again found that the ” critical chemical potential” is higher than the reasonable largest 
mass scale A of the NJL model thus the condition is nonphysical. In addition, it also de¬ 
mands a quite large and unrealistic electrical chemical potential Me- 

A comparison between the results with and without EN condition is shown in Table 3. 

Table 3 Comparison between the results with and without EN condition. 



Mci (MeV) 

(MeV) 

Mc 2 (MeV) 

EN 

345.6 

130.74 

656.87 

non-EN 

339.22 

102.37 

636.97 


It is seen that, in the case with EN condition, the first order critical quark chemical 
potential McI and the smaller least value point of V{m, fj,, fie) at the critical point are 
bigger than the corresponding ones in the case without EN condition. The former result 
is consistent with the conclusion derived in Ref. [5] that in the EN condition, a first order 
phase transition must happen at a larger value of Mu or m = M« + (2/3)Me than the one in 
the non-EN condition. In addition, with EN, the ’’second order critical chemical potential” 
Mc 2 to achieve m — >■ Mq is higher than the one without EN, thus it even more exceeds the 
largest mass scale A of the model. 

In short, inclusion of the EN condition does not change the qualitative behavior of chiral 
phase transitions of the model with the current quark mass Mq. Not only after a first order 
phase transition but also in the whole variations of to, the chiral condensates {qq) are always 
not equal to zeros within the frame of the NJL model where m < A. Therefore, we can not 
have complete restoration of the dynamical chiral symmetry breaking and finally, instead 
of talking about a second order phase transition, only get a crossover behavior of the order 
parameter to as m increases. 


4 A GENERAL ANALYSIS FOR ANY T AND ii CASE 


In this section, we will make a general analysis of the chiral phase transitions of the model 
in any T and m case. For simplicity and without loss of generality, we will consider only 
the case without EN condition. The effective potential V (to, T, m) at finite T and m in the 
mean field approximation may be expressed by [9] 


v(„,r,rt = - 12 


4Gs 


+Tln 




1 -I- I , E = + to^.(16) 


( 2 ^) 


Tin 


The gap equation dV{m, T, ^)/dm = 0 becomes 


TO 


Mq + 


12Gs 


TO 




1 


1 

^(E+fi)/T _j_ ^ 


(17) 


10 




















The second derivative of V(m, T, /i) over m may be expressed by 


d^V 

dm? 


1 

6 



P 


^ e(E-ti)IT I 
^{E-i£)IT 


[e(E-t^)/T + 1 ] 



( 18 ) 


By means of Eqs.(17) and (18) and the parameters given by Eqs.( 6 ) and (7), we may find out 
the least value points of V{m,T, fj,). For a given low T, analogous high /r chiral transitions 
to the ones at zero T as described in Sects. [2] and [3] will also appear and the discussions 
may be conducted similarly. Therefore, in Table 4 we list only the least value points of 
V (m, T, fi) at the critical chemical potential /r if a first order phase transition could happen 
at a given low T, or otherwise, the only minimal (least) value point of V (to, T, /i) at a higher 
T. 


Table 4 Variations of the least value points of V (to, T, /i) as change of T and fj,. 

For a given T and /i, toi and TO 3 represent the critical least value points with V (toi) = V{m 3 ) 
and TO 2 is the corresponding maximal value point. {V2)i and Vi (i = 1,2,3) represent re¬ 
spectively the values of d^V/dm? and V at the points mi. 


T 

(MeV) 

0 

20 

35 

37.1 

37.5 

50 

100 

100 

P 

(MeV) 

339.22 

335.52 

328.48 

327.27 

327.03 

327 

327 

500 

mi 

(MeV) 

301.47 

277.89 

219.93 

191.35 

186.23 

95.47 

44.81 

11.17 

(V2)i 

(MeV 2 ) 

Vi/lO® 

(MeV^) 

6920.41 

4419.43 

493.8 

7.287 

> 0 

> 0 

> 0 

> 0 

-6.915 

-7.100 

-7.5222 

-7.605 





m 2 

(MeV) 

211.71 

203.17 

189.32 

187.58 





iV 2)2 

(MeV 2 ) 

E 2 /IO® 

(MeV^) 

-2332.48 

-1591.39 

-220.27 

-4.605 





-6.858 

-7.074 

-7.5217 

-7.605 





m 3 

(MeV) 

102.37 

114.53 

156.79 

182.63 





(V^2)3 

(MeV 2 ) 

V 3 /IO® 

(MeV*) 

3042.59 

2211.86 

387.72 

9.19 





-6.915 

-7.100 

-7.5222 

-7.605 






It has been verified that at a given sufficiently low temperature T, as /r increases, we 
will first have a first order phase transition and at the first order critical quark chemical 
potential /r, we will always have the least value points of V{m,T, mi > m 3 ^ Mq, as 
shown in the left-first columns of Table 4. This implies again that after the first order phase 
transitions we have not achieved the condensates (qq) — 0 and the complete restoration of 
the dynamical chiral symmetry breaking, similar to the case of T = 0 . 

On the other hand, it is noted that as rising of T, for instance, from 0 to 37.1 MeV, 
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besides that the first order critical quark chemical potential /i will slightly go down (from 
339.22 MeV to 327.27 MeV), at the first order critical quark chemical potential /r, the 
bigger minimal point mi will go down (from 301.47 MeV down to 191.35 MeV) and the 
smaller minimal point m 3 will go up (from 102.37 MeV up to 182.63 MeV). This shows 
that the interval between mi and m 3 will become smaller and smaller. In the meantime, 
the curvature (oc \d'^V/drn?\) at the minimal points mi and m 3 (and also including at the 
maximal point m 2 ) will also become smaller and smaller. The final result will be that the 
two minimal points mi and m 3 will merge into a single one so that V (m, T, will have 
the mere minimal point mi left. This indicates that the first order phase transitions will 
end at r = 37.1 MeV and /r = 327.27 MeV. Once T > 37.1 MeV, the effective potential 
V (m, T, /i) will always have the only minimal (also least value) point left whose location mi 
will smoothly decrease as a further increase of T and fi, as shown in Table 4. As a result, 
in the T — fjL phase diagram of the model, T = 37.1 MeV and /i = 327.27 MeV will become 
a critical end pint i.e. at this point the first order phase transition line ends. 

It should be emphasized that for the appearance of the critical end point, a decisive factor 
is the nonzero-ness of the chiral condensates {qq) i.e. m 3 Mq after a first order phase 
transition. It is the nonzero-ness of {qq) that makes it possible that the two minimal points 
mi and m 3 could be moved and finally merge into a single one. Otherwise, if it is assumed 
that we may have m 3 —>■ the fixed Mg i.e. {qq) —>■ 0 after a first order phase transition, and 
furthermore, the only minimal point mi of V (m, T, ^) at a higher T can finally be reduced 
to Mg (similar to a variation of the order parameter in a second order phase transition), 
then we will obtain, similar to the chiral limit case with Mg = 0, a tricritical point, rather 
than a critical end point. 

Of course, if we have m 3 ^ Mg after a first order phase transition but assume that as 
a further increase of T and/or both m 3 and the least value point mi of V{m,T, ^) at a 
higher T could be reduced to Mg, then besides that the explicit chiral symmetry breaking 
is always kept in the phase transitions and the critical end point remains to exist, we will 
also get a ’’second order” phase transition line. However, we will prove that this situation 
can not practically happen within the frame of the NJL model. 

As is shown in Table 4, for the least value points m of V(m,T, ^), we always have 
m = m 3 > Mg after a first order phase transition or m = mi > Mg at a higher T. In 
addition, similar to that indicated in the right-last columns of Table 4, m will decrease as 
a further increase of T and /r. Thus a natural question is that if at some high value of T 
and /i, the limit m —>■ Mg, corresponding to total disappearance of the chiral condensates 
or equivalently, complete restoration of the dynamical chiral symmetry breaking, could be 
achieved within the frame of the NJL model. For giving the answer of this question, we 
must solve the gap equation (17) with m replaced by Mg. In this case Eq.(17) becomes 



^(Eo-n)/T _|_ 2^ g{Eo+^i)/T _|_ 2 



= 0 , 


(19) 


which contains only the parameters A and Mg, independent of G 5 . The problem is reduced 
to find out the solutions of Eq.(19) about T and /r. 

First for a given T, consider possible values of /i. 

We begin with T = 0. When T —?> 0, we have -|- 1] —>■ 0, thus Eq.(19) will 

be reduced to Eq.(9) which has the solution fj, = as given in SectHJ 

For a finite T, we have numerically solved Eq.(I9) and obtained the results listed in 
Table 5. 
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Table 5 Solutions of Eq.(19) for some finite T. 


T (MeV) 

0 

20 

50 

100 

200 

p (MeV) 

636.97 

1102.02 

1858.53 

3138.57 

5781.7 


These data indicate that at these (and believably any) finite T, the resulting values of /i 
are always bigger or much bigger than A, thus are not physically reasonable. 

On the other hand, we may consider the solutions of Eq.(19) when /i < + Mg. 

In fact, numerical calculations give T = 4.06 x 10® MeV for all /r < a/A^ + Mg. In view 
of the working precision of the calculations, this merely indicates that T —^ oo is the only 
solution of Eq.(19) if /r < A^ + Mq . However, the result T —>■ oo is obviously physically 
unrealistic. 

The above results can be given a simple physical explanation. It is noted that the gap 
equation (17) at finite T and fi is identical to 


TO = Mo - 2Gs{qq)T, 


( 20 ) 


where {qq)T is the thermal quark-antiquark condensates which can be obtained from the 
quark-antiquark condensates (qq) in vacuum through replacement of the quark propagator 
by its form in thermal field theory [la¬ 
in fact, the quark-antiquark condensates (qq) in vacuum can be expressed by 


{qq) = NfNc 


d'^p i 

__ 

(27r)‘* ^ — TO -I- ie ’ 


( 21 ) 


where Nf and Nc are the quark’s flavor and color number respectively. Taking the real-time 
formalism of thermal field theory, then in Eq.(21) making the replacement 


> — m + ie 


COS^ 9r) 


I — m + ie 


-—h sin^ 9„ 


siM 9p = 0(/)n(/ -pi)+ 0(-/)n(-/ + p) 


— m — le 
0 


- d) = 


1 


e(p°-G/T + 1 ’ 

we may express the thermal quark-antiquark condensates by 


{qq)T = --iNfNcm 


d'^p 

Wr 

d^p m 




:-27rsin^ 9p5{p^ — to^) 


= -ANfNc J [1 - n{E - p)- n{E + p)] 


6 to 


E 


dp 


1 - 


1 


1 




( 22 ) 


(23) 


where Nf = 2 and Nc = 3 have been taken. In view of Eq.(23), we immediately see that 
Eq.(20) is just the gap equation (17). 

Now it is easy to find the difference between Mq = 0 and Mq ^ 0 for Eq.(20). When 
Mq = 0, i.e. in the chiral limit, Eq.(20) will become 


.. 2Gs 
1 =- {qq)T- 

TO 

It may have the solution to = 0 at some finite and reasonable T and p which is identical 
to that the thermal condensates (gg)T = 0 and indicates a complete restoration of the 
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dynamical chiral symmetry breaking. On the other hand, when Mq ^ 0, if one also expects 
such complete restoration of the dynamical chiral symmetry breaking, then it must have the 
result that 

= 0 - 

However the proceeding discussions in this section just show that when the current quark 
mass Mo ^ 0 the above result could not be attained at a finite T and a reasonable ^ < A 
i.e. within the frame of a NJL model where p,, T < A. We indicate that it is just a nonzero 
{qq)T that can have a continuative smooth decrease as rising of T and ^ thus could lead 
to the crossover behavior of the order parameter m indicating chiral symmetry breaking at 
high T and high /r. Otherwise, if m —?> Mq i.e. the condensates (qq) = 0 could be achieved 
finally within the frame of the NJL model, then, as stated above, it will lead to a ’’second 
order phase transition line” rather than the crossover behavior of m in the T — fj, phase 
diagram. 

In brief, it is the nonzero-ness of the chiral condensates induced by the current quark 
mass Mq which leads to the appearance of the known critical end point and the crossover 
behavior of the order parameter m in the phase diagram of the NJL model. In view of the 
simulating feature of the NJL model for QCD, such phase diagram is very well qualitatively 
consistent with the ones obtained based on real QCD dynamics with the current quark mass 

M- 

In addition, the condensates {qq) 7 ^: 0 at a finite T, ^ < A demonstrated here will also 
force one to consider the interplay between the quark-antiquark condensates {qq) and the 
diquark condensates {qq) when one researches the color superconductors in zero T and 
middle fj, case, where fj, has exceeded the first order critical quark chemical potential fj,ci 
and the diquark condensates could be formed. This point has been noticed and discussed 
earlier [Hllli], and the analysis made in the present chapter certainly strengthen theoretical 
grounds of the above consideration. 

5 CONCLUSION 

In this chapter, by means of a quantitative investigation of the order parameters in the 
ground states we have reexamined and analyzed systematically the effect of the current 
quark mass Mq on nonzero-ness of the chiral condensates {qq) in chiral phase transitions 
of a two-flavor NJL model simulating QCD respectively in the cases of zero T and high /i 
without and with the electrical neutrality condition as well as in the case of any T and /i 
without the electrical neutrality condition. The nonzero-ness of the chiral condensates {qq) 
induced by the current quark mass Mq is found to have a double meanings: one is that the 
condensates {qq) keep to have quite large values after a first order phase transition at a low 
T and a large /r = /ici; the other one is that although at a high T or a large fi > /ici, the 
chiral condensates {qq) will smoothly decrease as a further rising of T and/or /i, they can 
never be reduced to zeros within the frame of the NJL model, i.e. in the condition when 
/i, T < A, where the momentum cutoff A of the loop integrals must be viewed the largest 
physical mass scale of the NJL model. 

For the latter one, in fact, the mathematical solving of the gap equation indicates that 
in the case of T = 0, the order parameter m —>■ Mq i.e. the chiral condensates {qq) = 0 
could be achieved only at /r = a/ A^ -|- Mq if the EN condition is not imposed, and even at a 
higher ^ if the EN condition is imposed. In the general case, it has been proven that when 
the current quark mass Mq ^ 0, the limit m —>■ Mq can be attained only in the conditions 
that ^ > OY ^ A^ + Mq if T = 0 and finite, or T —00 if /r < A^ + Mq. Obviously, 
these conditions are all non-physical and unrealistic. 

The above results indicate that once the current quark mass Mq exists, it is impossible 
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to achieve complete restoration of the dynamical chiral symmetry breaking within the frame 
of the NJL model. This includes both after a first order phase transition and in the whole 
process of chiral phase transitions. Theoretically this reflects a close and deep interconnec¬ 
tion between the explicit and the dynamical (spontaneous) chiral symmetry breaking. 

The nonzero-ness of the chiral condensates {qq) plays a decisive role in the structure 
of the phase diagram of the NJL model with the current quark mass Mq. It is just the 
nonzero-ness that leads to the appearance of a critical end point rather than a tricritical 
point, and a crossover behavior of the order parameter at high T and /i rather than a second 
order phase transition line in the known phase diagram of the NJL model. In addition, it 
also gives a strong theoretical grounds for the point of view that one must consider the 
interplay between the quark-antiquark condensates {qq) and the diquark condensates {qq) 
when researching the color superconductor at zero T and middle fi. 

To sum up, the analysis made in present chapter of the nonzero-ness of the chiral con¬ 
densates induced by current quark mass and its physical effects may certainly deepen our 
theoretical understanding of the mechanism of chiral phase transitions of the NJL model. 
It is also an excellent example from which one can understand how a source term of the 
Lagrangian can greatly affect dynamical behavior of a physical system. 
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